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Abstract 

We consider the computational complexity of Hamiltonians which 
are sums of commuting terms acting on plaquettes in a square lattice 
of qubits, and we show that deciding whether the ground state mini- 
mizes the energy of each local term individually is in the complexity 
class NP. That is, if the ground states has this property, this can be 
proven using a classical certificate which can be efficiently verified on a 
classical computer. Different to previous results on commuting Hamil- 
tonians, our certificate proves the existence of such a state without 
giving instructions on how to prepare it. 



1 Introduction 

Understanding the ground state properties of spin systems on a lattice is of 
central importance in many-body physics, but at the same time, it is a highly 
challenging problem in many scenarios. An important step in understanding 
its difficulty has been the insight that computing e.g. the ground state energy 
of a classical spin system is, in general, an NP-complete problem [1]: While 
the energy of any given spin configuration can be easily computed, finding 
the configuration with minimal energy is in general a difficult task - it can 
be as hard as any problem in NP, i.e., any problem whose solution can 
be efficiently verified. For quantum spin systems, an additional difficulty 
arises: Generally, we cannot even expect to have an efficient description 
of the ground state. Thus, it seems that the only statement we can make 
about the difficulty of the problem is that given a quantum register with 
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the ground state, we will be able to efficiently estimate its energy using 
a quantum computer. Indeed, it has been shown that this is the best we 
can say, as the problem of estimating the ground state energy of a quantum 
system is a complete problem for the class QMA ^tSj, the quantum analogue 
of NP: It contains all problems which have a quantum solution which can 
be efficiently checked on a quantum computer, and thus, determining the 
ground state energy of a quantum spin system is as hard as any of these 
problems; in fact, the problem retains its hardness even when restricted to 
two-dimensional lattices of qubits with nearest-neighbor interactions [1] or 
one-dimensional chains [5]. 

It is an interesting question to understand the reasons underlying the 
additional complexity of quantum spin systems as compared to classical 
systems. To this end, restricted version of the problem which lie between 
classical and general quantum spin Hamiltonians have been studied: For 
instance, it has been shown that so-called stoquastic Hamiltonians, where 
all off-diagonal elements are negative, have a complexity which lies in be- 
tween NP and QMA, as those systems can be related to classical random 
processes [SI [7]; in fact, these are exactly the Hamiltonians which allow for 
Quantum Monte Carlo simulations as they do not exhibit a sign problem. 

Another restricted class of Hamiltonians are commuting Hamiltonians, 
that is, Hamiltonians which can be written as a sum of mutually commuting 
few-body terms. For those systems, all terms can be simultaneously diago- 
nalized, just as for classical systems; however, the corresponding eigenbasis 
can be highly entangled, making it unclear whether a useful classical descrip- 
tion of the ground state can be provided. In fact, commuting Hamiltonians 
encompass systems which exhibit rich non-classical physics, in particular 
models with topological order and even anyonic excitations, such as Ki- 
taev's toric code and quantum double models [8], or Levin and Wen's string 
net models [9] . Commuting Hamiltonians are also of interest since the fixed 
points of renormalization flows in gapped phases are expected to be commut- 
ing Hamiltonians, and thus understanding their structure might give insight 
into the structure of gapped quantum phases. Finally, understanding the 
complexity of commuting Hamiltonians is of interest in quantum complex- 
ity, as it might be a step towards a quantum PCP theorem, which would 
assess how the difficulty of estimating the ground state energy is related to 
the desired accuracy which is integer for commuting projectors. 

What is know about the complexity of finding the ground state energy 
of commuting Hamiltonians, or rather, of determining whether the ground 
state minimizes each term in the Hamiltonian individually - the commut- 
ing HAMILTONIAN problem? For lattices in two and higher dimensions. 
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COMMUTING HAMILTONIAN is an NP-hard problem, as it e.g. encompasses 
the Ising model on a planar graph [1]. On the other hand, it is not clear 
whether the general commuting hamiltonian problem is inside NP, since 
it is not clear in general how to provide an efficiently checkable description 
of the ground state. For two-local (i.e., two-body) Hamiltonians, Bravyi 
and Vyalyi [TO] have shown that the problem is in NP by using C*-algebraic 
techniques (their result also implies that one-dimensional commuting Hamil- 
tonians are efficiently solvable); subsequently, Aharonov and Eldar [llj have 
proven containment in NP for Hamiltonians with three-body interactions 
both for qubits on arbitrary graphs, and qutrits on nearly-Euclidean inter- 
action graphs. In both of these works, the classical certificates do not only 
prove the problem to be in NP, but can in fact be used to construct constant 
depth quantum circuits which prepare the ground state. This, in particular, 
implies that the corresponding Hamiltonians - including qutrits with three- 
body interactions - cannot exhibit topological order I12j. On the other 
hand, Kitaev's toric code, which is the ground state of a commuting Hamil- 
tonian with four-body interactions of qubits, does have topological order, 
and thus, we cannot expect any approach which yields a low-depth circuit 
to work beyond three qutrits. 

In this paper, we study the commuting hamiltonian problem on a 
square lattice of qubits with plaquette-wise interactions, and prove that it 
is in NP. That is, we consider a square lattice of qubits, with a Hamilto- 
nian with mutually commuting terms acting on the four qubits adjacent to 
each plaquette, and show that the problem of deciding whether its ground 
state minimizes the energy of each local term in the Hamiltonian is in NP: 
i.e., in case the ground state has this property, a classical certificate exists 
which can be checked efficiently by a quantum computer. Our approach dif- 
fers considerably from the aforementioned approaches in that the certificate 
cannot be used to devise a quantum circuit for preparing the ground state, 
and is thus also applicable to systems with topological order; it should be 
noted that the same holds true for the proof in Ref. [TO] that commuting 
HAMILTONIAN with factorizing projectors is in NP. 

2 The setup 

We will consider a 2D square lattice with spins on the vertices, and either 
open or periodic boundary conditions. The Hamiltonian 

H = Y,h, 
p 
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consists of plaquette terms hp which act on the four spins adjacent to the 
plaquette p, and we impose that ah the terms in the Hamiltonian commute, 
[hp,hq] = \/p,q. 

As the Hamihonian terms commute, there is a basis of eigenstates of H 
which are simultaneously eigenstates of all the hp. We would like to know 
the computational difficulty of the following problem, called commuting 
HAMILTONIAN: Is there an eigenstate {ip) of H which minimizes the energies 
for all hp individually, i.e., are the ground states of H also ground states 
of each hp? In the following, we will show that in the case of qubits, the 
existence of such a state can be proven within NP, i.e., there is a classical 
certificate which proves the existence of such a \^), and which can be checked 
efficiently classically. Note that on the other hand, it is clear that the 
problem is NP-hard, as it e.g. encompasses classical Ising spin glasses in 
a field which are known to be NP-hard even for two-level systems [1]. 

For the following, it wih be useful to reformulate commuting hamilto- 
nian as follows: Define the local ground state projectors Hp as the projectors 
onto the ground state subspace of hp] the Hp commute again, [11^,11^] = 0. 
Then, 

ncs = n 

is the projector onto the subspace spanned by the states which are ground 
states of ah hp. Since commuting hamiltonian asks whether such states 
exist, it is equivalent to asking whether IIgs 0. 

3 Commuting Hamiltonian in NP 

3.1 Two layers 

We start by coloring the plaquettes of the square lattice black and white in 
a checkerboard pattern, and denote the set of black and white plaquettes by 
Vb and Vw, respectively. Let 

Ub= Yl^p and Uw= n 

P&Vb P&Vw 

be the projectors onto the joint ground state space of the black and white hp, 
respectively; then, commuting hamiltonian corresponds to determining 
whether IIbIIvi/ ^ 0. This is equivalent to asking whether 

iiiHB^w] / (1) 
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(this can be seen using eigendecompositions of Ilw ^-nd Hb), and we will 
consider this formulation of the problem from now on; to prove COMMUTING 
HAMILTONIAN is contained in NP, we therefore need to show that a classical 
certificate for the validity of ([T]) can be provided. 

A helpful example to keep in mind in the first part of our discussion is 
Kitaev's toric code [8]: There, Up = i(l + Z®^) for p & Vb, and Up = 
i(l + X®^) for p G Vw, with X and Z the Pauli matrices. 

3.2 The structure of one layer 

In the following, let us study the structure of each layer individually (we 
will w.l.o.g. choose black). To this end, we will use a result of Bravyi and 
Vyalyi based on C*-algebraic techniques [10]; a detailed explanation of those 
techniques can also be found in [11]. The basic insight from Ref. [10] is 
the structure of two commuting projectors. Consider two projectors L = 
Lab ® Ic and i? = 1^ ® Rbc acting on a space Ha ^T-Lb ® H-c', the two 
operators overlap on Hb- Now consider the Schmidt decompositions 

L = Y,^L^Bi and R = Y,B'r0C},, 

i.e. tr[j4^(A-^)'''] = for i / j, and similarly for B\^, and C]j. Then, 
[L,R] = implies that [B\^,B-'^] = for all i,j, and thus, and B]^ 
span commuting C*-algebras, cf. [10]. Using the standard form of finite 
C*-algebras, it follows that the space Tis has a canonical decomposition 

nB = ^n®n%^n''z (2) 

— ■Hg 

where the B\^ span the full matrix algebra on Ti'^ while acting trivially on 
the rest, and correspondingly for B^j^ and V.'^. 

This shows that the space Hb can be cut into direct sum "slices" (the 
a) such that in each slice, L and R act on independent subsystems. More 
formally, there exists a decomposition 1 = X^^vr^ oH-Lb, with tTq the pro- 
jectors onto 11%, such that 

[tTq,, L] = and [tTq, R\ =Q 

and thus 

L = ^ TTaLnp = ^ TTaLTTa , R = ^ TTaRlTa , 
a,/3 a =-ia o =Jla 
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where L'^ and i?" act on different subsystems Tif^ and H'^, i.e., factorize. 
Note tliat tfie above decomposition allows to compute the vTq, and thus the 
La and Ra efficiently. 

Each vertex in the black sublattice is acted upon by exactly two com- 
muting projectors 11^; thus, we can apply the preceding argument to all 
vertices to find decompositions vr^^, tt^^ = 1, of the Hilbert space at 
each vertex v, such that projected onto the slice a = {av)v^v factorizes, 

(g)n^^=n<„nB<„ 

(this implicitly defines the Hp), and 11^ can be written as 

Hb = (g) n^^ = 5; (g) (3) 

Note that in the sum on the right hand side, we implicitly regard the tensor 
products as being canonically embedded into the full Hilbert space, and we 
will use this convention in the following. 

The analogous decomposition can be performed for the white sublattice, 
yielding a (in general different!) decomposition 

n^ = © (g) ® 4- (4) 

p p£Vw P&Vw 

Note that in order to distinguish the decomposition for the the black and 
the white layer, we strictly use labels a and /3, respectively; moreover, we 
denote the projectors decomposing the white layer by vf^^. 

E.g., for Kitaev's toric code the vr^^ are projectors onto the Z eigenstates, 
and the vf^^ onto the X eigenstates. 

3.3 Combining the layers 

Using the stucture of He and Hw, Eqs. ^ and (jl]), we can rewrite the 
COMMUTING HAMILTONIAN problem, Eq. ([T]), as 



0/ J^tr 

d,i3 



(2) n?)((8) n^p 

P&Vb P&Vw 



(5) 



recall that we regard the tensor products as being canonically embedded 
into the full Hilbert space. Since each of the traces is positive (as it is the 
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trace of the product of two positive operators), the above is equivalent to 
the existence of a and /3 such that r2(a, /3) ^ 0, where 



n{d,(3) :=tr 



(2)n?)((8) n^^ 

P&Vb P&Vw 



(6) 



Thus, we can ask the prover to provide us as a certificate with some a 
and (3 for which rj(a, /3) ^ 0; if we can moreover show that rt{a,(3) can 
be computed efficiently (or rather in NP), this will prove that commuting 

HAMILTONIAN is in NP. 

Note that ^2(5, /3) can be interpreted as the overlap of the (unnormal- 

ized) states (^11^ and (^IIp, both of which are tensor products of states 
supported on individual plaquettes, but with different partitions in the two 
layers. Computing such an overlap can in general be as hard as contracting 
Projected Entangled Pair States (PEPS) [l3], i.e., PP-hard [H]: any PEPS 
can be expressed this way by using one layer for the bonds and the other 
for the PEPS projections. Of course, the fact that these states arise from 
two commuting layers 11^ and liw yields additional constraints, and we will 
show in the following that those constraints allow for the efficient evaluation 
of il(a, /3) in the case of qubits. 

For Kitaev's toric code, e.g., we could choose a = /3 = (0, . . . , 0): This 
yields 

n{aj) = tr[(|0)(0|)^^(|+)(+|)n = 2-^ / , 

proving the existence of a zero-energy ground state; note that this certificate 
does not carry any information on how to prepare the state. 



3.4 Computing the overlap 

Let us now show that for a lattice of qubits, the overlap 0(a, /3) can be 
computed efhciently. To this end, we will show that the computation of the 
overlap can be decomposed into a product of overlaps of one-dimensional 
structures which can be computed efficiently. 

Let us first consider the black layer. For each vertex v, the decompo- 
sition ([2]) of the local Hilbert space can either be trivial (no direct sum) 
or non-trivial, 1 = X] • the former case, this implies that at most 
one of the adjacent plaquette terms Yip acts non-trivially of vertex v; in 
the latter case, the sum consists of exactly two one- dimensional projectors 
TT^^, making use of the fact that the Hilbert space at each site is a qubit, 
i.e., two-dimensional. We will denote the set of vertices with a non-trivial 
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decomposition in the black layer by J-b, and in the white layer by J-y/ (in 
which the one-dimensional projectors are labelled vf^^). 

In 0(a, /3), all vertices in Tb U contribute only a one-dimensional 
subspace and thus can be traced out: That is, all vertices in Fb n Fw 
can be removed (taking care whether the overlap of the one-dimensional 
projectors is non- vanishing) , while for vertices where only one layer has a 
one-dimensional decomposition, this yields new effective plaquette terms pp 
in the other layer by projecting the original plaquette terms lip onto that 
one-dimensional subspace; thus, the problem of checking whether il(a, /3) 
is non-zero reduces to computing the overlap of the new effective plaquette 
terms pp. Formally, this reads 



rj(a,/3) = tr 



((2)n? 



peVw 



(A) 



tr 



( n n 



TTr- 



{*) 



(B) 



v£Tb<~^J^w 



((g) Pp){<S) Pp 

peVw 



I peVB 

Here, we have used in step (A) that for all v G Fw, 



(7) 



(8) ^^P=^l{<S) 



P&Vw 



and anologously for the black layer; in step (B), we have defined new effective 
plaquette terms pp by virtue of 



Pp — trjTguJ'w 



P£Ve 



n 

weTw\J^B 



TTf- 



(2)n? n n 

P&Vb w&Tw\^b 



and correspondingly for the white plaquettes; the first factor in ([7|) takes care 
of the terms in Fb n Fw in (*)• Note that the pp are now only supported 
on those vertices not in Fb U Fw-, as those have been traced out. 

The task of checking whether n{a,j3) ^ has thus been reduced to 
checking this for ([7]): For the first term, this can be clearly done efficiently. 
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Figure 1: Computing the overlap Eq. ^ for qubits. The diamonds 
at the vertices of the square lattice denote the qubits. The connected 
black dots mark qubits on which the pp act non-trivially on plaquette 
p. If the Pp of two plaquettes act non-trivially on the same qubit, we 
say that they "overlap" ; note that this cannot happen for diagonally 
adjacent plaquettes. Overlapping ppS form structures which we need 
to contract to compute the overlap Eq. Q . a) Patterns forming one- 
dimensional chains can be contracted efficiently, as the size of the 
boundary stays constant for any contiguous region, b) Branching 
structures do in general not allow for efficient contraction. However, 
we show that such structures cannot occur, by proving that for any 
plaquette p, pp can overlap non-trivially at most with two adjacent 
Pp/'s. 



and so, it remains to prove that the overlap 



can be computed efficiently. Note that the pp are now supported on pla- 
quettes of a square lattice with vertices missing. Moreover, while the pp do 
no longer commute, in each layer at most one pp acts non-trivially on each 
vertex; we will make use of this fact repeatedly in the following. 

The situation encountered in computing the overlap Q is depicted in 
Fig. [TJ Here, the dots in each plaquette denote the vertices on which pp 
acts non-trivially (the lines just connect the vertices involved in pp). If the 
Pp on adjacent plaquettes act non-trivially on the same qubit (we will say 
they "overlap" ) , they form connected structures which we need to contract in 
order to evaluate 0. For one-dimensional structures as the one on in Fig. [1^, 
this contraction can be carried out efficiently: One starts from one pla- 
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quette and proceeds along one direction of the one-dimensional chain, always 
tracing out the degrees of freedom on the inside. This way, at every point in 
the computation only the state at the boundary (which has fixed size) needs 
to be stored, and thus, the contraction can be carried out efficiently. On the 
other hand, branching structures like the one in Fig.[Tl3 can in general not be 
contracted efficiently, since the size of the boundary is a priori not bounded; 
in fact, e.g. quantum circuits, and even postselected quantum circuits, can 
be encoded this way, making such contractions in general a computationally 
hard task. 

However, as we will show in the following, the structures formed by the pp 
in 0, Eq. ([8]), will always be one-dimensional, and thus Q can be computed 
efficiently. To this end, we will consider the state pc on a plaquette C (the 
"central" plaquette), and show that it can overlap non-trivially with the 
states Pp of at most two of the adjacent plaquettes, thus ruling out branching 
structures as the one on the right of Fig. [TJ We will make intensive use of the 
fact that in each layer, at most one plaquette term pp can act non-trivially 
on any given vertex; in the graphical notation of Fig. [H we will highlight 
this fact by placing a cross opposite of any dot: 



Here, the dot indicates that the corresponding pp acts non-trivially on a 
vertex qubit (the diamond in the center), while the cross indicates that the 
corresponding pp does act trivially. Note that this in particular implies that 
pc can at most overlap non-trivially with the four horizontally and vertically 
adjacent plaquettes from the other layer, but not with diagonally adjacent 
plaquettes. 

The simplest case is when the state pc on the central plaquette C involves 
only two vertices non-trivially, for instance 



Now, both qubits 1 and 3 can be acted upon non-trivially by at most one 
white plaquette - the other has to be empty, following the rule ([9]) that 
opposite of any dot there has to be a cross; this way, only one-dimensional 




(9) 
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structures can arise: 




This clearly holds for any possible pc which acts non-trivially on only two 
qubits, and for any configuration of the adjacent plaquettes; it follows that 
only one-dimensional structures can emerge this way. 

In order to understand the cases where pc acts non-trivially on three or 
four qubits, let us first analyze the following situation: 




Here, pc acts non-trivially at least on qubits 1 and 2, and pL acts non- 
trivially on qubit 1 which implies that pr acts trivially on it. In the following, 
wc will show that this implies that pT also has to act trivially on qubit 2. 
We will prove this by contradiction, so assume that pT acts non-trivially 
on qubit 2. Since pT is obtained from the original projector 11^ on that 
plaquette by a partial projection on some of the other vertices, this implies 
that acts non-trivially on qubit 2 (where it spans the whole C*-algebra, 
since we have traced out all vertices where this was not the case). On the 
other hand, pi and thus Hl acts non-trivially on qubit 1, and thus spans 
the whole C*-algebra on it. Since 11^: and IIt commute, this means that 
I\.T acts trivially on qubit 1; that is, IIt and He need to commute on qubit 
2 alone: However, since Ht spans the whole C*-algebra on qubit 2, this 
would imply that He and thus pc has to act trivially on qubit 2, giving a 
contradiction. Thus, we have the following "Lemma" : 




Let us now consider the case where the state on the central plaquette 
involves all four qubits non-trivially, and let us start by assuming w.l.o.g. 
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that PL acts non-trivially on qubit 1: 




This imphes that qubit 1 is not acted upon by pT- Using Eq. (jlOp . we infer 
that pt cannot involve qubit 2 either, 




Eq. (jlOp also shows that ps has to act trivially on qubit 4 - otherwise, pi 
would act trivially on qubit 4 and thus qubit 1, which it doesn't: 




In order to obtain a branching structure, both ps and pR need to act non- 
trivially on some of the qubits. However, if ps acts non-trivially on qubit 
3, Eq. (jlOp implies that pR has to act trivially on qubits 2 and 3: 




This shows that for pc acting non-trivially on all four qubits, the central 
plaquette can only couple to at most two adjacent plaquettes, forming one- 
dimensional structures. 
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It remains to study what happens in the case of tripartite entanglement 
on the central plaquette. We start from the following configuration: 




(Note that we don't make any assumptions on how the Hilbert space of qubit 
2 decomposes.) Clearly, in order to obtain a branching structure, either pT 
has to act non-trivially on qubit 1, or pn has to act non-trivially on qubit 
3. We consider w.l.o.g. the first possibility and infer from Eq. pop that pi 
has to act trivially on both qubits 1 and 4: 




In order to obtain a branching structure, we now have to get both pn and 
Pb involved. However, if pn acts non-trivially on qubit 3, Eq. (fTOl) implies 
that Pb has to act trivially on both qubits 3 and 4, and thus, the structure 
formed around pc will again be one-dimensional. 

Together, this shows that the overlap 0, Eq. (I8|), decays into one-dimen- 
sional structures for which the overlap can be computed efficiently. In turn, 
this implies that for given a and /3, VL{a, j3) can be computed efficiently, and 
thus, the commuting Hamiltonian problem on a square lattice of qubits with 
plaquette interactions is in NP. 



3.5 Finite accuracy 

In the preceding proof, we have assumed infinite accuracy, but as we will 
now show, our argument still applies if we compute with finite accuracy. To 
this end, let denote the number of qubits in the system; we will need 
to show that the computation time scales as poly(iV). We assume that the 
Hamiltonian terms are given exactly and can be represented with poly(A^) 
digits. First, note that the trace in Eq. ([T]), which equals the sum in Eq. ([5]), 
evaluates to an integer, and thus, there exists at least one pair (a, /3) such 
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that Q{a,f3) > 2~^^. If we request this particular (a, /?) as a certificate, it 
is sufficient if we can determine r2(a, /3) to 2N + 1 = poly(A^) digits. (This 
is crucial, since there can be (a, /3) for which 0(a, /3) is arbitrarily small.) 
0(a, /3) is obtained from contracting a polynomial number of terms which 
are either Hp (which are known exactly) or tt^^ and 7f^^, and the latter can 
be determined to poly(A^) accuracy from the C*-decomposition ([2]), which 
is the solution to a (fixed-size) eigenvalue problem. It follows that il(a, /?) 
can be computed to the required poly(A^) accuracy in poly(A^) time, and 
our proof still applies. 

4 Summary 

We have studied the commuting hamiltonian problem on a square lattice 
of qubits with plaquette-wise interaction and shown that the problem is N P- 
complete. Differently speaking, we have shown that there exists a classical 
certificate for the fact that the ground state of the system minimizes each 
term locally which can be checked efficiently on a classical computer. The 
central idea for our proof has been to split the system into two layers in each 
of which the commuting terms overlap on individual sites, and to argue 
that the existence of a state minimizing all local terms is equivalent to 
the existence of a pair of ground states for the two layers with non-zero 
overlap. Each layer could be decomposed using the C*-algebraic techniques 
introduced to the problem in [lOj, allowing to find an efficient description 
of its ground state subspace. Finally, we showed that the overlap of ground 
states of two layers can be computed efficiently, by showing that it gives 
rise to of one-dimensional structures only. A somewhat surprising feature 
of our approach is that while it certifies the existence of a ground state, 
it cannot (to our knowledge) be used to devise a way how to prepare the 
ground state; in fact, due to the possibility of having topological order in 
such systems, any circuit preparing their ground states would need to have 
at least logarithmic depth, or linear depth if it was local |12j . 

Our method does, in principle, also apply beyond qubits: We can still 
split the system into two layers, decompose both of them into direct sum 
slices Q and /3, and ask the prover to provide labels a and /3 with non-zero 
overlap i}{d,f3). While we cannot make sure any more that 0(a, /?) can be 
computed efficiently, we can always ask the prover to also provide us with 
an instruction on how to efficiently contract the states, in case there is a 
way to do so, e.g. by providing the optimal contraction order. In particular, 
this applies to the case where the decomposition in the direct sum gives 
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one-dimensional spaces, such as in Kitaev's toric code or quantum double 
models; as well as to cases where the pp are separable states. Our idea also 
applies to any other graph which can be split into two layers in such a way 
that the C*-technique of [lOj can be applied to each of them, and in fact to 
any type of Hamiltonian which is composed of two layers with eigenbases for 
the zero-energy subspace whose overlap can be computed efficiently, such as 
for product bases. Note that on the other hand, a decomposition into three 
layers cannot be used for our purposes, since for three positive operators A, 
B, and C, tr [^i3C] can have both real and imaginary parts of either sign, 
so that Eq. ([5]) is no longer equivalent to Eq. ([6|) being non-zero (as the Hp 

and Hp do not commute any more). 

An interesting open question relating to the present approach to the 
problem is whether it can be generalized beyond qubits. For four-level sys- 
tems and beyond, this is likely not the case, since the local Hilbert space can 
decompose into two qubits, and thus operators commuting on a single spin 
can both act non-trivially on it, i.e., Eq. ([9|) does not hold any more. On 
the other hand, for qutrits this is not the case once we have fixed a slice in 
the direct sum; yet, it is not clear how to establish a version of Eq. (llOp . In 
particular, the non-trivial projections vr^^ can now have both rank 1 and 2, 
and in the latter case we cannot simply trace out the corresponding degree 
of freedom; it is however not clear that this does rule out an analogue to 
Eq. (do]). 
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